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Abstract. 

Kilohertz Quasi-Periodic Oscillations (QPOs) have been detected in many accreting X-ray binaries. It has been suggested 
that the highest QPO frequencies observed in the modulation of the X-ray flux reflect a non-linear resonance between two 
modes of accreting disk oscillation. This hypothesis implies certain very general predictions, several of which have been 
borne out by observations. Some of these follow from properties of non-linear oscillators, while the others are specific to 
oscillations of fluid in strong gravity. A 3:2 resonant ratio of frequencies can be clearly recognized in the black-hole as well 
as in the neutron-star QPO data. 
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1. INTRODUCTION 

The resonance model of Kluzniak and Abramowicz 1 1 ] explains pairs of high frequency QPOs observed in neutron 
star and black hole (microquasar) sources in Galactic low mass X-ray binaries as being caused by a 3:2 non-linear 
resonance between two global modes of oscillations in accretion flow in strong gravity 1 . In this report, we shortly 
review rather convincing observational arguments that support the resonance model: 

1 . The rational frequency ratio 3:2. For microquasars , the observed ratio v u / V] of the upper and lower frequency is 
3/2. For neutron star sources the ratio varies in time, but its statistical distribution peaks up, within a few percent, 
at the 3/2 value. 

2. The inverse mass scaling. In compact sources, v u ~ 1 /M. 

3. The frequency -frequency correlation. The upper and lower frequencies of a particular neutron star source are 
linearly correlated along the "Bursa line", v u —Av\+B, with A ^ 3/2. Observational points occupy a finite 
sector of the Bursa line which typically crosses the reference line v u = (3/2) Vi at the "resonance point" . 2 

4. The slope-shift anticorrelation. For a sample of several neutron star sources, the coefficients A ,B of the Bursa 
lines corresponding to individual sources in the sample are anticorrelated, A = 3/2 — Z?/(600Hz± Av), where 
Av <C 600 Hz is a small scatter. 

5. Behavior of the rms amplitudes across the resonance point. For the few sources examined so far, the difference 
between rms amplitudes of the lower and upper QPO changes sign at the resonance point. 

One should have in mind that except the inverse mass scaling the above properties are typical of a non-linear 
resonance in any system of two weakly coupled oscillatory modes. They give no direct clues on how to answer several 
fundamental questions connected to the specific physics of the oscillations of accretion disks: 

• Which are the two modes in resonance? An often discussed possibility is that they are the epicyclic modes [see, 
e.g-,|5l], but there are some difficulties with this interpretation. 



1 For a recent review of observations see 0] and for a review of the theory see 01. 

2 Only a few examples of sources for which the actually known data do not cross the 3:2 line are known [see, e.g.,0|. 



• What are the eigenfrequencies of the modes in question ? It was recently found that pressure and other non- 
geodesic effects have stronger influence on the frequencies of some modes than previously thought |6, 7]. If this 
question is answered, one may be able to precisely determine characteristics of strong gravity, in particular the 
black hole spin (see Section l2~2l . 

• How are the modes excited? What is the energy source that feeds the resonance? It was suggested that the 
excitation in the neutron star case could be by a direct influence of the neutron star spin (e.g. due to a magnetic 
coupling) |8, 9], and in the black hole case by influence of the turbulence 111 (ill 1 111 211 . Understanding of this point 
may provide strong constrains on the magnitude and nature of turbulence in accretion disks. 

• How is the X-ray flux modulated by accretion disk oscillations? It was suggested that in the case of neutron star 
sources modulation occurs at the boundary layer lfl3ll . and in the black hole sources it is due to gravitational 
lensing and Doppler effect — at least partially [see, e. ^..[l4lfT^I . The mechanism of modulation must be related 
to connection between QPOs and the "spectral states" lUal . 



2. BLACK HOLE HIGH FREQUENCY QPOS 
2.1. The 3:2 ratio and 1 /M scaling 

The 3:2 ratio is a strong argument in favor of explaining the QPOs in terms of a resonance [see, e.g.,yi for details]. 
In general relativistic phenomena around a mass M, all frequencies scale as 1 /M. This agrees with black holes having 
lower QPO frequencies than neutron stars and is suggestive of a relativistic origin for kHz QPOs [1]. McClintock and 
Remillard |25] found that in microquasars the frequencies do indeed scale inversely with mass: 

V " =2 - 793 (S ^ (1) 

The 1/M scaling could be used to estimate the black hole mass in AGNs and ULXs [;27|, if HF QPOs were to be 
discovered in those sources (see Figures^^. 



2.2. Measuring the black hole spin 

According to the epicyclic resonance hypothesis, the two modes in resonance have eigenfrequencies v ra d (equal to 
the radial epicyclic frequency) and v v (equal to the vertical orbital frequency Va or to the Keplerian frequency Vk). 
Several resonances of this kind are possible, and have been discussed [see, e.g.,|5|]). 

Formulae for the Keplerian and epicyclic frequencies v ve i-t and v ra( j in the gravitational field of a rotating Kerr black 
hole with the mass M and internal angular momentum a (spin a) are well known |29J: 

v 2 = vl(l -6x~ l +8ax- 3/2 -3a 2 x- 2 ^ , 

V 2 e = v| (l -4ax- 3/2 + 3a 2 X - 2 ^ , 



r r 
X ~~ix;~ 2GM/c 2 1 



TABLE 1. Frequencies of twin peak QPOs in microquasars. 



Source* 


v u [Hz] 


AVupp [Hz] 


V,[Hz] 


Avi[Hz] 


2v u /3v! - 1 


Mass 1 " [M ] 


GRO 1655^0 


450 


±3 


300 


±5 


0.000 


6.0 — 6.6 


XTE 1550-564 


276 


±3 


184 


±5 


0.000 


8.4 — 10.8 


H 1743-322 


240 


±3 


166 


±8 


-0.036 


not measured 


GRS 1915+105 


168 


±3 


113 


±5 


0.009 


10.0 — 18.0 



* Twin peak QPOs first reported by ITllTllT9t 120.121 . 

f See 1 22, 23, 24, 25] for the microquasars. Note that there is also the different estimate for GRO 1655^40: 
M= (5.4±0.3)M7El. 
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FIGURE 1. The microquasar kHz QPOs frequencies in a 3:2 ratio 1 25] may be extrapolated to the black hole in the Galactic 
center I2HI2I. 
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FIGURE 2. Left: 1/M scaling in black hole sources. The solid line is the best 1/M fit for microquasars by McClintock and 
Remillard 1 25] and possibly identified with the 3:2 epicyclic resonance frequency for a ~ 0.97, dotted line is plotted for the same 
resonance and a = 0. Right: Position of the center of a torus fulfilling the 3:2 resonance condition, as a function of the the /3 
parameter charecterizing non-slenderity of the torus In this Pseudo-Newtonian result, the resonant frequencies increase with /3. 
A similar effect in the Kerr metric would lower the black hole spin parameter estimated for a given QPO frequency. 



For a particular resonance n:m, the equation 

nv r = mv v ; v v = Vq or v K (3) 

determines the dimensionless resonance radius x n - m as a function of spin a. Thus, from the observed frequencies and 
from the estimated mass one can calculate the relevant spin of a central black hole 113(113 ill . 

Results of this procedure were summarized, e.g., in 1331 . Several resonance models give the values of spin in the 
range a 6 (0, 1). In particular, the 3:2 epicyclic parametric (or internal) resonance model supposed to be most natural 
in the Einstein gravity 11311 implies microquasar black hole spin a ~ 0.9. However, the most recent results of the spin 
estimate for GRO 1655-40 fe.g. . 13311 claims the spin to be ~ 0.7. It could be interesting that recently proposed 3:2 
periastron precession resonance 1 34] implies the spin of GRS 1915+105 to be also a ~ 0.7. 

On the other hand, above mentioned resonance spin estimates are based on a resonance between epicyclic frequen- 
cies equal to those for exactly geodesic motion Recently |6, 7], it was found that pressure and other non-geodesic 
effects have stronger influence on the frequencies than previously thought. At the moment, it is difficult to give plausi- 
ble estimate how this effects change the spin predictions as the theory is not developed enough, but work on this is in 
progress. A rough estimate, based on Pseudo-Newtonian calculations of ||7j] and the assumption of the Kerr case being 
qualitatively same, says that corrected values of spin required from 3:2 parametric resonance model will be rather 
lower then was supposed (see Fig. [2] right panel). 



3. NEUTRON STAR HIGH FREQUENCY QPOS 



As was shortly reminded above, for the black holes v u and Vi appear to be fairly fixed having the well defined ratio 
V u /V| = 3/2. Contrary to this, V u , Vi in neutron star sources vary by hundreds of Hertz, along "Bursa lines" 

v u =Avi+B. (4) 




V] [Hz] 



FIGURE 3. Bursa plot: The linear fit of the neutron-stars data (dashed line) is obviously inconsistent with the 3:2 relation, but 
the frequency ratio is peaked to the 3:2 value, i.e. close to the 3:2 resonance 1 35]. Psaltis et al. 1 36] has shown, in a different context, 
that frequency pairs of all sources cluster along a si ngle line. Bursa (2002, unpublished) also pointed that each source follows a 
slightly different line. It was demonstrated by 1 37, 38, 3 9] that variations of v u and V] along the line can be explained within the 
non-linear resonance model for QPOs (see Section l3~2l . 



3.1. Ratio vs. frequency distribution 

The distribution of the V u /vi frequency ratio in neutron star sources was shown by Abramowicz et al. 1 35] to cluster 
near the 3:2 value. This result was later criticised by Belloni et al. [40] who thought that the relevant quantity must be 
the distribution of a single frequency. 

We stress that this is a misunderstanding and these two distributions touch different problems. Linear correlation 
between the frequencies doesn't imply the ratio clustering because for this the location of data segment is important 
as well (see Figures|4]and[5}- Nevertheless, Belloni et al. 14011 confirmed the 3:2 clustering in the neutron star sources 
and shown also the very interesting clustering close to the related resonant frequency. 




FIGURE 4. A simulated segment of 57 frequency-frequency points (related by v u = 0.75vi +450) with a flat distribution in the 
lower frequency is placed approximately at the position of Sco X-l data in v-v diagram. Resulting histogram of ratio is peaked at 
the 3:2 value. 




3.2. Anticorrelation between slope and shift 

Abramowicz et al. 114 ill have shown that linear relations @ for twelve neutron star sources are anticorrelated. Below 
we shortly remind the connection between this anticorrelation and properties of weakly coupled nonlinear oscillators 
[see^, for details]. 

The frequencies of non-linear oscillations may be written in the form 

Vi = V,° + Avi , Avi = V,° ( K X a\ + K u al) , 

v u = v u ° + Av u , Av u = v°(X l af + X a a 2 tt ), (5) 

where Vj and Vy are the eigenfrequencies of oscillator and a/, a u are the amplitudes of oscillations. 

If the two amplitudes, a\ and a u are correlated, i.e., if they are functions of a single parameter s, one gets a linear 
relationship between v u and Vi, obtained from lowest order expansion 3 in the parameter s. 

The inferred coefficients A and B of linear relation @ will vary from source to source, but if the eigenfrequencies of 
oscilations (Vj , V„) are universal for the sources (being fixed by the common space-time metric and mass), then these 
linear relations themselves will satisfy 

v° 1 1 

A = ^-^ B=A °-^ B - (6) 

Therefore, the slope A and the shift B are anticorrelated. If the eigenfrequencies of oscillators slightly differ of factor 
<7 (e.g., due to the difference of the mass and connected 1 jM frequency scaling) but the ratio of eigenfrequencies is 
fixed, this anticorrelation reads 

A=A Q --^B. (7) 

Thus, for given type of resonance, the individual pairs A, B should locate on lines inside a triangle in the A-B plane 
with a quite well determined vertex at [0, Aq], and with the size of its base proportional to the scatter in Vi. In particular 
for the 3:2 resonance, the vertex should be [0, 1 .5]. Figure [6] suggests that this is indeed the case. 

In principle, for the 3:2 parametric epicyclic resonance model the eigenfrequencies of oscillations can be directly 
identified with the epicyclic frequencies of Keplerian motion at the location of this 3:2 resonance. This implies that 
the generic neutron star mass connected to the straight line in Fig. [6] is about 1M (first noticed by Bursa 2003, 
unpublished). However, recent results 0] suggest that due to a pressure and other non-geodesic effects, this mass will 
be higher. 



3 Note that if one includes more terms in the expansion, the frequency-frequency correlation could deviate from straight line. 
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FIGURE 6. The coefficients of the linear frequency-frequency relation for kHz QPOs in twelve neutron-star sources: 

1) 4U 0614, 2) 4U 1728, 3) 4U 1820, 4) 4U 1608, 5) 4U 1636, 6) 4U 1735 EJ; 7) 4U 1915 EJ; 8) XTE J1807 |4]; 9) GX 17+2 |44]; 

10) GX 34+0 1 45]; 11) GX 5-1 0; 12) ScoX-1 HJ]. 

The data are consistent with a linear relationship intercepting the A axis at the value 1.5 (Note good quality of fit by straight line). 
The scatter (shaded area) expressed in terms of a coefficient from eq. is most likely from the interval a £ (0.8, 1.2). Under the 
assumption that difference in the eigenfrequencies is caused by different neutron star mass, the mass should differ of factor ~1.5. 



3.3. RMS amlitude evolution across the resonance point 

Quite recently, Torok discovered the strong change in character of the observed QPO amplitudes (Torok et al. j47ll ). 
The difference between the rms amplitude of the lower QPO and the rms amplitude of the relevant upper QPO, i.e., 
the quantity 

AA,-„„(v) = A v r ] m {v) -A^(v), where v = Vi (8) 

is rather well correlated with frequency and changes its sign across the 3:2 resonance line which is illustrated in 
Figures 0and[8] 

It was noticed by Horak that such phenomenon resembles general behavior of mechanical systems crossing an 
internal resonance (Horak et al. |48]). As some frequencies approach a rational ratio the modes involved in the 
resonance start to exchange energy. The slow motion through the resonance perturbs the strict balance between both 
directions of energy exchange and finally causes increase of one amplitude at the expense of the other one. 
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FIGURE 7. The difference of the rms amplitudes vs. lower QPO frequency for two neutron star atoll sources 4U 1728-34 and 4U 
1608-52. Each datapoint corresponds to one continuous segment of observation (error bars without central circle denote insignificant 
observations). Cyan points correspond to the difference in double peaks. Orange points follow from single QPOs identified through 
the Q-factor [see 49], their counterparts are here supposed (probably not quite realistically) to have rms equal zero. It seems that 
up to the 'resonance point' the amplitudes are larger in the case of the upper frequency, equal when 'Bursa line' pass 3:2 ratio and 
above this, the lower QPO has the stronger amplitude. Based on mAl in preparation. 
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FIGURE 8. Left: The difference (yellow curve and symbols as in previous figure) and total (grey curve) of the rms amplitudes for 
atoll source 4U 1636-53. The continuous curves are calculated from polynomial interpolation of the values of rms 1 50] for the upper 
and lower frequency obtained through shift and add over all PDS. Right: Double peaks vs. single QPO which are distinguished 
using Q-f actor and placed in accord to linear fit. Up to the intersection of the fit with the 3:2 line, lower frequencies are often 
supressed (i.e., difficult to find) as the single 'upper' frequency dominate, while above this the situation is opposite. Based on [47] 
in preparation. 
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